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ABSTRACT 



(N 

I . Context. We present a detailed study of the X-ray properties of the galaxy cluster Abell 548b (z=0.04), using XMM- 

O ' Newton data, and discuss the connection between the thermal properties and the presence of two extended relic radio 

■ ^ ' sources located at the cluster periphery. 

C/^ , Aims. We wish to analyze the dynamical state of the cluster and confirm the presence of a major merger. We will discuss 

C^ . the merger effects on the extended nonthermal emission. 

Methods. We analyzed the X-ray brightness distribution, the surface brightness profiles in different directions, and the 
spectral properties in several cluster regions. Moreover, to better understand the dynamical history of this cluster, we 
performed an optical analysis of the cluster galaxies. 

Results. From the analysis of the temperature distribution and of the surface brightness profiles, we find evidence 
of a shock in the northern cluster region, just before the location of the two extended peripheral relics. From the 
optical analysis, we find that the cluster galaxies show a large av, and two components are needed to fit their velocity 
distribution. Observational results were compared with a cluster simulation. The maps of gas temperature and density 
distribution from the simulation agree with the observational data in the case of a cluster merger nearly perpendicular 
^^ ■ to the plane of the sky and in the state after the maximum core collapse. 

^^ ' Conclusions. We conclude that we are observing a galaxy cluster in a major merger phase, just after the maximum core 

OO ! collapse. The mass ratio is about 1:2, and the merger collision is nearly perpendicular to the plane of the sky. A shock is 

^^ , present in the northern cluster region, and it is very likely responsible for the electron reacceleration and the magnetic 

field amplification that will give cause the cluster relics. The relative position of the shock and the relics is strongly 
affected by projection effects. 
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1. Introduction spect to the cluster center, and another possible relic source 

of small size (C), ne ar in projection to the cluster center 



Cluster formation is a still active process in the local (iFeretti et al. l[2QQ6h . These diffuse radio sources have no 

Universe. Galaxy clusters exhibit small-and large-scale sub- obvious connection to cluster galaxies, but are rather asso- 

structures that clearly indicate of merging activity with ciated with the ICM 
other clusters or groups of galaxies. The X-ray observations 
of the thermal gas in galaxy clusters allow us to study the 

hot intra cluster medium (ICM), which is the main baryon Diffuse cluster radio sources appear to fall into two cate- 

source in the Universe. Thanks to the XMM-Newton and gories: halos, which are centrally located in the cluster, rel- 

Chandra observations, we can obtain detailed temperature atively regular in shape, and unpolarized, and relics which 

distributions of the ICM and understand the influence of are peripherally located, fairly elongated, irregular and of- 

the cluster's formation history on the thermodynamics of ten highly polarized. The presence of these large regions 

the IGM. (300-600kpc), which show a diffuse synchrotron emission, 

The galaxy cluster Abell 548b (z =0.04) is a merging reveal the existence of a large-scale magnetic field and of 

cluster with two extended radio relics (A and B), located relativistic particles in the cluster volume. These nonther- 

in a cluster peripheral region, on the same side with re- mal components have earned a lot of attention in the last 

years since their formation and evolution is connected to 

Send offprint requests to: L. Solovyeva, the cluster's dynamical history, and therefore they are im- 

e-mail: lilia.solovyeva@cea.fr portant for a complete description of the ICM. 
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Fig. 1. X-ray image from XMM data in the 0.3 - 4.5 keV 
energy band, from EPIC cameras. Green contours are the 
radio emission at 1.4 GHz; A, B, and C indicate the radio 
reUcs. 



Present models of the origin of relics suggest that 
they are tracers of merger shocks. Relativistic particles 
are accelerated by shocks either by diffusive shock accel- 
eration or adiabatic recomp ression of fossil radio plasma 
(|Ensslin and Bruggen I l2QQ2:Hoeft & Bruggen 2007). The 
production of outgoing shock waves at the cluster periph- 
ery is indeed o bserved in nume rical simulations of cluster 
merger events (|Schindler 11200 2*). 

The cluster A548b is located in the cluster system A548, 
which includes several substructures: A548a, A548b, and 
A548 -2 associated with optical condensations (|Davis et al. I 
Il995f ). The combination of X-ray and optical data show 
that A548 is a cluster that is still in a formation phase and 
that has a complex morpho logy. T he ROSAT data, pre- 
sented by iDavis et al. I (|l995l ) and by iNeumann fc Arnaud I 
(|1999[ ), indicate that this cluster is in a merging phase, 
hence not yet dynamically relaxe d. A cooling core is clearly 
not present f lWhite et allll997[ ). 

In this paper we present a study of the dynamical state 
of the galaxy cluster A548b by detailed X-ray analysis 
with new XMM- Newton data to understand the connec- 
tion between X-ray properties and the presence of radio 
relics. In Sect. 2 we present a summary of the radio prop- 
erties. In Sect. 3 we present the new data, their reduction, 
and the results of the detailed image and spectral analy- 
sis (surface brightness profiles, 2D /3-model analysis using 
XMM-Newton and ROSAT data, total temperature pro- 
file, temperature profiles in different regions, temperature 
map). Section 4 gives the relationship between X-ray and 
radio data. In Sect. 5 we compare our results with optical 
data and simulations to analyze the dynamical status of the 
cluster. Results are discussed in Sect. 6, and conclusions are 
found in Sect. 7. 

We assume a concordance cosmology with Hq = 
70kms-^Mpc-\ Qrn = 0.3, Qa = 0.7. At the redshift of 
the Abell 548b cluster (z = 0.04), the angular scale of V 
corresponds to a linear size of 47 kpc. 




Fig. 2. X-ray corrected image from XMM data with chosen 
sectors, smoothed with a Gaussian with FWHM = 25^'. 
Superimposed are sectors discussed in Sect. 3.3.3 

2. Summary of radio properties 

Figure [T] shows the contours of the radio emi s sion a t 1.4 
GHz, obtained and discussed bv I Feretti et al. I (|2006l ). We 
concentrate here on the extended diffuse radio sources. A 
diffuse source, named A, was detected in the western clus- 
ter region. This extended source has quite a regular mor- 
phology and is not associated with any galaxy. Also diffuse 
emission, named relic B, was detected in the northern clus- 
ter region, around a strong point-like radio source. This 
extended emission is filamentary and irregular. The em- 
bedded radio galaxy is not at the center of the diffuse radio 
emission, but located at its southern boundary. The two 
diffuse radio sources A and B are strikingly similar in the 
total flux density at 1.4 GHz and with similar size. The 
relics A, B are polarized at a level of 30 per cent. No in- 
formation about Faraday rotation is available. They have 
steep spectra, with a spectral index value of ~-2± 1. 

Another apparently diffuse source (labeled C) was de- 
tected closer to the cluster center, to the northwest, at 
about V from the galaxy triplet VV 162, classified in the lit- 
erature as a dumbbell galaxy and coincident with the X-ray 
centroid. The diffuse radio source C is elongated with no ob- 
vious optical identification, and it is polarized at a 7% level. 
At high resolution, it shows a central diffuse emission and 
two brighter and more compact spots possibly identified 
with cluster galaxies. No sign of a radio source coincident 
with the brightest cluster galaxy is detected. The total flux 
density of this complex in the high-resolution image is only 
about 1/3 of that in the image at lower resolution, thus in- 
dicating a diffuse structure of lower brightness. Therefore, 
this source has been classified by Feretti et al. (2006) as a 
relic. 

3. X-ray analysis 

3.1. Observations and data reduction 

Abell 548b was observed on 2006 February 17 with 70 ksec 
for EPIC MOS cameras and 64 ksec for pn camera. Figure 
[1] shows the A548b observation of XMM-Newton with the 
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radio relics contours (A,B,C). The pointing position was 
chosen to be between the cluster center and the position 
of radio relics to have information on the cluster periph- 
eral region and to study the connection between the X-ray 
and radio emission. For t he data analysis we used the back- 
ground of J. Nevalainen (|Nevalainen et al71l2QQ5h. We also 
checked our results using the background of iLumb et al. I 
(|2QQ 2) and obtained similar results for temperature and 
surface brightness profiles. We used the niethod of double 
background substraction bv lArnaud et al. I (|2Q02l ) (see next 
section). For the analysis, we used the XMM- Newton data 
from EPIC cameras (MOS 1,2, and pn) and the XMM- 
Newton Science Analysis System (SAS) for data reduction. 
In the MOS 1,2 data set we took into account event pat- 
terns to 12 and in pn data - patterns to 4, flag=0. In 
the spectral analysis, we used the data from three EPIC 
cameras, while we used only the EPIC cameras M0S2 and 
pn in the image analysis and 2D beta analysis, because we 
did not have any statistics in CCD6 for camera MOSl. For 
background observation, the sky coordinates in the event 
files were modified using the aspect solution of the clus- 
ter observation. From the data counts, we detected and 
excluded the periods of "flares". The flnal observing time 
after flare subtraction was 45 and 33 ksec respectively. We 
normalized the background with the count rate of observa- 
tion in the high energy bands and exposure time of obser- 
vation. For correction of the vignetting effect in observa- 
tion and background dat a, we used the weighted function 
(|Maierowicz et al.l |2QQ2[ ) . The extended cluster emission 
was clearly detected, and we analyzed the data in the 0.3- 
4.5 keV energy band to optimize the signal-to- noise ratio. 
Visible point-like sources were masked with circles, and ex- 
cluded from the observational data in our spatial and spec- 
tral analysis. We filled these gaps with the mean value of 
surrounding environment by multiplying by random values 
from Poisson distribution. 



3.2. Double background subtraction 

After cleaning the flare events, the XMM background is 
dominated by the cosmic X-ray background (CXB) and non 
X-ray background (NXB) induced by high-energy particles. 
The background subtraction is very crucial when reducing 
extended diffuse sources. To prope rly remove it, we use d 
the double-subtraction method by lArnaud et al. I (|2QQ2I) . 
It consists of two steps: (i) we subtracted the correspond- 
ing normalized blank fleld obtained using the same spa- 
tial and energy selection (NXB component), and (ii) then 
we subtracted the residual components using the data in 
the outer part, outside the cluster (CXB component). The 
CXB component of the background is affected by vignetting 
of the X-optics. In t he flrs t step we used the background 
of iNevalainen et aLl (|2QQ5l ). In the second step, we mini- 
mized the x^ ill the outer regions. The background counts 
were normalized using the count rate of our observation in 
the high-energy bands and taking the observation exposure 
time into account. 



3.3. X-ray brightness distribution analysis 

3.3.1. X-ray image 

In Fig. [1] we show the A548b XMM- Newton image (colors) 
uncorrected for background or exposure. The superposed 



contour represents the radio e mission. The radio re lics are 
named A, B, and C according to 'Ferett i et al. I (12006'). From 
this image, it is clear that the cluster emission it extend out 
to 4' - 5' from the central region in north direction (N) and 
up to 3' - 4' in south direction (S). The pointing position 
(see Sect. 3.1) is far from the X-ray peak in the XMM image, 
and the cluster center is strongly affected by the vignetting 
effect. For this reason we used the position from the ROSAT 
image ([Bohringer et al. Il2004i) as the cluster center. 



3.3.2. Corrected X-ray image with modeling of background 
from observation data 

In the data analysis of extended X-ray emission from clus- 
ters of galaxies, the background estimate and subtraction is 
very important. In our case we also have the problem that 
the vignetting effect can seriously affect our data because 
of the peripheral pointing center. 

For all the 2D analysis, we needed to obtain a 
vignetting-corrected and background-subtracted image. 
This is specially important since A548b was observed in the 
border of the fleld of view where vigneting effects is more 
signiflcant. Since one component of the background, the 
CXB, is vignetted by the X-rays optics when the other, the 
NXB, is not, we need to estimate their relative contribution 
to be able to subtract them correctly. In the outer region, 
where the cluster contribution could be neglected, we can 
model the observation by BKG = CXB * 1/F -h NXB. 
Since the vignetting function, VF, is known from calibra- 
tion, one can disentangle the CXB from the NXB contri- 
bution by fltting this model to the background data in this 
external region. From this model and our outer region pro- 
flle, we minimize the y^ to achieve the two components 
feolovvevaet al. II2007D . 

The obtained image of background was subtracted, and 
we flnally corrected our image for the gap and the pres- 
ence of point sources. We fllled the regions where the mask 
value is zero with the mean value of the surrounding envi- 
ronment by multiplying by random values from Poisson dis- 
tribution. We then summed the images for cameras M0S2 
and PN, excluding the MOSl camera because of the ab- 
sence of CCD6. Figure [2] shows the cluster image after the 
background subtraction, with vignetting, gap, and point 
sources correction. This image also shows the chosen direc- 
tion of sectors for study, and the image is smoothed with a 
Gaussian fllter of a = 7 " . 



3.3.3. Surface brightness profile 



Abell 548b is a non relaxed cluster ([Davis et al. 1119951 ), so 
we did not expect it to show a gas distribution consistent 
with a single /3-model. We thus obtained flts with a /5-model 
in different directions. We performed a PSF convolution on 
the beta- model in the flt of the surface brightness proflle. In 
the obtained image for each XMM camera (using standard 
data reduction, with the background of J. Nevalainen), we 
summed and rebinned the photons into concentric annuli 
with a size of ~ 3 .3 ^^ centered on t he ce nter of the X-ray 
emission given bv iBohringer et alTI (|2004f ) (see sect. 3.3.1) 
for each camera. In Fig. [3^1 we show the total radial surface 
brightness proflle and the obtained /3- model flt. 

We also derived the surface brightness proflles' averag- 
ing data in three sectors: N, Western (W), and Eastern 
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Fig. 3. a) Total surface brightness profile and /?- model fit. b) surface brightness profile in a different direction from 
the center cluster N, W, E and comparison with /3-model fit obtained from radial surface brightness profile, c) surface 
brightness profile only in north direction from center cluster in particular: NE, NW and between relics A and B and 
comparison with /?- model fit obtained from radial surface brightness profile. 



(E), defined as in Fig. O The results are in Fig. [8)3. As 
a second step to obtain more detailed information in the 
relic's direction, we obtained surface brightness profiles in 
sectors NE, NW, and in between the two relics (FiglSh). 
The different sectors, except those between the two relics, 
are shown in Fig. [2l (Note the sector between relics A and 
B is the angle 45-70 degree if the sector N is the angle 10- 
110 deffle e.) W e fitted each surface brightness profile with a 
/3- model (|Cavaliere and Fusco-Femiano l il976). We did not 
use the first ten points to avoid contamination from the 
galaxy triplet VV 162, present at the cluster cluster. We 
note that the brightest cluster galaxy is not identified with 
any discrete radio source and that no optical counterpart 
is associated with the extended emission identified as relic 
C (|Ferettiet al. II2006I and Sect. 2). 

The results of the /3- model fit are given in Table [TJ All 
the fits at any position angle are poor, as indicated by the 
high values for the reduced x^. This agrees with the fact 
that A548b is a non relaxed cluster. We als o compared our 
results with those of iNeumann fc Arnaud I ([1999) obtained 
from ROSAT, and the comparison confirms that the region 
where data show the more significant perturbation is at ~ 
1.5' and at ~ 5', in each of the directions for our research. 



Table 1. A548b, /?- model fit in different directions 



direction 



c(arcmin) x /d.o.f 



total 
West 
East 
North 



0.67: 



7+U.U8 
-0.02 

0.56+°-°^ 



3.01] 



0.80 
0.62 



0.03 
+ 0.06 

0.07 
+ 0.03 

0.03 



North-West 
North-East 
Between relics A B 



Neumann& Arnaud 
from 1' to 5' 



P, 70+U.U6 

'-'•'^-0.07 

60+°°^ 

U.OU_Q Q3 

71 7+0-10 



1+0.13 
-0.20 
9 ^9+0.44 

OQ-,+0.40 

9 74+0.20 



1+0.42 
-0.47 

^+0.22 
-0.25 

O-O0_0.47 



171/40 
96/33 

78/34 
134/40 



3.79: 

2.22^ 



161/39 

65/38 
50/32 



0.52zb0.06 2.82±0.04 



54+0-02 



9 91+0.14 

^•^-'--0.14 



1.35 
43/30 




Fig. 4. Residual image obtained from corrected XMM- 
Newton image, after subtraction of the 2D P model. 

Table 2. Best-fitting results of the 2D /?- model 



Parameter 


Best-fitting 






XMM 


Rosat 


RciC) 


2.46 


2.76 


Rc2n 


2.87 


3.34 





0.69 


0.65 


Pa 


1.92 


1.99 


RA 


05:45:29 


05:45:28 


Dec 


-25:55:56 


-25:55:58 



3.3.4. 2D-;^ analysis, X-ray residuals 

To better analyze the brightness distribution in the cluster, 
we performed a two-dimensional fit to the cluster brightness 
using a modified 2D- /3-model that allows for two different 
core radii along the two principal axes of the cluster im- 
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age ellipse. We estimated the deviations from the model, in 
the inner 6'. Figure [4] shows the residuals after the model 
subraction from the XMM- Newton Gaussian-smoothed im- 
age. The residual map is smoothed with a Gaussian of 25 " 
FWHM. If one computes the significance of the structures 
of this size, we found roughly 10 a for each structure inside 
a 10 arcmin distance from the cluster center. 

The same analysis was performed in the ROSAT data 
and we obtained similar results. In particular from both 
data sets we obtained a very similar residual map show- 
ing a clear decrease in the emission at 1.2', and 7', and 
an increasing emission at 4.5' around the cluster center, in 
the relic direction (N and W). This result agrees with the 
surface brightness profiles obtained at the different direc- 
tions (Sect. 3.3.2). The best fitting parameters are listed 
in Table [2l Note that XMM data give better results then 
ROSAT data since the XMM image was obtained after a 
more careful subtraction of the background (CXB and NXB 
components). 

The present profile analysis confirms that A548b is a 
perturbated cluster, with a complex dynamical history. 
Moreover, we observe an asymmetrical perturbation from 
the 2D- /3- model displaced from the cluster center. To bet- 
ter understand the dynamical state of A548b, we need to 
perform a spectral analysis now. 



3.4. Spectral analysis 

To perform the spectral analysis of the XMM-Newton 
data, we used three cameras of XMM-Newton, the back- 
ground from J. Nevalain en, and the method of double back- 
ground substraction by lArnaud et al. I (|2002') , and we fit- 
ted the spectrum with XSPEC using a redshifted MEKAL 
lasma emission model with absorbtion Nh = 2 • l{fi^cm~^ 



pla^ 



Davis et al. 1995) and abundance = 0.3. Reported errors 



in the text and tables are at a 68% confidence level. 

The best fit to the mean cluster temperature in the inner 
3' gives a central cluster temperature kT = 3.4 ±0.1 keV 
and a reduced x^ = 1.02. 

We checked our results using models with free abun- 
dance, and obtained similar results and, in particular, same 
temperature values. 



3.4.1. Temperature in different regions 

To study the influence of the ICM gas on the radio relics, 
we extracted the X-ray spectrum in different regions near 
the radio relics. Chosen regions are shown in Fig. [5l The 
results are shown in the Table [3l A temperature higher 
than the cluster average derived in the previous subsection 
was found in regions 6 and 4, located between the cluster 
center and relics A and B, with temperatures of 4. 6^0.55 ^^^ 
4.5^0 3 keV, respectively. The fits in the regions 1 and 2, i.e. 
in the relic regions A and B, are not conclusive because of 
the poor statistics in this peripheral cluster area. However, 
in region 1 (relic A) we obtained a temperature of 6.21^9 
keV, suggesting a possible temperature increase in the relic 
location. To study the influence of the ICM on relic C, we 
extracted the spectra in the regions 7 and 3. The obtained 
temperatures are 3.010.2^^^ ^^^ ^-^-o.s keV, respectively, 
so very close to the cluster average temperature in the inner 
3' ( 3.4l^;l keV). 




Fig. 5. Chosen regions of A548b for spectrum extraction 



Table 3. Spectral fits results by regions 



region 



nkeY)/xted 



1 

2 
3 

4 
5 
6 

7 


6.23l?-^/L06 
3.74lJ-^/0.89 
3.95t°-^/0.80 
4.46t°-^/0.88 
4.04l°-2/0.88 
4.62l°-^/0.93 
3.00t°-^/0.91 


cluster 


2' 
3' 
5' 


3.28t°-}/0.93 
3.37+°-}/1.02 
3.59+°-}/1.09 




Fig. 6. a) Temperature map of A548b obtained by sec- 
tors with contours of relics, b) "Significance map" = 

J- ohs J- average 
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3.4.2. Temperature Map 

The temperature map is the best measurable indicator of 
equilibrium in a system. We obtained temperature maps 
by integrating by sectors or boxes {2' in size) and by an X- 
ray wavelet spectral mapping algorithm. In all three cases 
we obtained similar results and find a clear temperature 
increase at ~ 5' from the cluster center followed by a tem- 
perature decrease at ^ 7'. The most significant temperature 
increase was from 3.2keV to 6.5 keV, obtained in the relic 
direction on the 2-3 keV. Figure [6] we show the temperature 
map obtained in different sectors where the temperature 
trend is more visible. In Figl6)3 we have the 'significance 
map' that shows the difference between observed and aver- 
age temperature divided on temperature error, where aver- 
age temperature is 3.3 keV through the whole the cluster. 

3.4.3. Temperature profiles in sectors 





Fig. 7. a) temperature profiles of A548b in different direc- 
tions: N, S, E, W b) temperature profiles of A548b in relic 

directions: N, NW, NE. 



For a better understanding of the previous result, we 
derived temperature profiles in the different sectors for the 
directions: N, S, E, W (as in Fig. [2]). We extracted spectra 
in each sector up to 10' from the cluster center, and the 
size of sectors was chosen similar to the size annuli used to 
derive the radial temperature profile . 

Figure [7^ shows obtained temperature profiles in each 
direction. A significant increase in the temperature at about 
5^ for the N direction, (in the region just before relics A 
and B), is clear. By splitting the N sector in the two sub- 
sectors NW and NE, we observed the same increase of the 
temperature in the region between 4' - 6^ before the relic 
location (see Fig. [JId). The highest temperatures are about 
5-6 keV . 

3.5. Do we detect shocks toward the relics? 

Figure [8] shows the surface brightness and temperature pro- 
files in the direction of the radio relics A and B. We clearly 
observe an increase at 5' in the temperature profile (see Fig. 
Uh ) and a significative drop at 7' in surface brightness, just 
before the relic location. On one hand, if one considers that 
statistics are much better in brightness and that surface 
brightness is less affected by the averaging along the line 
of sight, one should use this profile to estimate the shock 
position and its mach number. On the other hand, the de- 
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Fig. 8. Obtained perturbation from X-ray data in the sur- 
face brightness profile and in the temperature on the N 
direction to relics A and B. 



tection of a clear jump in the temperature could also be 
interpreted as a shock signature. That the projected posi- 
tions of the shock are different in the two profiles should 
then be interpreted as different projection/averaging effect 
for surface brightness and temperature (also taking into ac- 
count that we need a larger bin to estimate the temperature 
and then suffer stronger effects). Due to these statistical 
considerations, it seems better to define the shock position 
using the surface brightness profile. We, thus compute the 
temperature profile again so as to have bins just before and 
after the shock (see Fig. [8]). 

Considering a possible shock in the N direction and 
taking the surface brightness profile into account (see Fig. 



^+0.33 
-0.27 



at 



[8|), the density compression given by — is 1.74 

7^ We can use this value to constrain the kinematics of 
the merger: the Mach number M of the shock can be de- 
rived by employing the Rankine-Hugoniot shock relations 
(iLandau fc Lifshitz I [19591 . iMarkevitch fc Vikhlinini I2QQ7D 
for density or temperature jumps of 



P2 _ (1+7)M^ 
pi 2+(7-l)M2 

T2 _ (27)M2 - (7 - 1))((7 - 1)M^ + 2) 
Ti (7 + 1)2M2 



(1) 



(2) 



where the subscript 2 denotes material behind the shock 
(the "postshock" region), and 1 denotes material ahead of 
the shock (the "preshock" region). For the intracluster gas, 
we used 7 = 5/3. 

Using the above compression value and Eq. 1, the data 
agree with the presence of a moderately supersonic outflow 
with a Mach number M= 1.5 ±0.1. This result implies a 
temperature jump of T2/T1 = 1.54, in very good agreement 
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with the measured discontinuity, at 1' : T2/T1 = 1.41q"4. 
In a first conclusion, keeping in mind both the different 
averaging and projection effect along the line of sight for 
temperature and density and the fact that both jumps lead 
to similar mach number, we think that we probably detect 
a relatively mild shock toward the relics. 



4. Relationship between X-ray and radio data 

4.1. X-ray versus radio emission 

Current models propose that merger shocks are at the ori- 
gin of relics. In our case, the two radio relics A and B are 
well outside the cluster core radius (^100-150 kpc), at pro- 
jected distances of 500 and 430 kpc from the cluster cen- 
ter, at the boundary of the detected X-ray emission. The 
poor statistics do not allow us to derive conclusive infor- 
mation on the temperature in the relic regions, whereas a 
significant temperature increase is detected in the region 
before the relics. The presence of perturbations in the X- 
ray brightness distribution and of variations in the cluster 
temperature is a clear indication of a disturbed cluster, thus 
supporting the connection between the formation of relics 
and a strong cluster dynamical activity. On the other hand, 
it seems that there is no point-to-point coincidence between 
the radio regions and the shocked region. We show that it 
is possible to explain the link between the radio emission 
and the X-ray properties by invoking projection effects. 

The small relic C, located at about 1' from the cluster 
center, i.e. within the cluster core radius, in the NW direc- 
tion, is not associated with any obvious X-Ray substructure 
or temperature jump. This source consists of three com- 
ponents: a central diffuse emission and two brighter and 
more compact spots possibly associated with cluster galax- 
ies (see Figs. 6 and 7 and the discussion in Feretti et al. 
2006). The total flux density present in the high-resolution 
images is only about 1/3 of the flux density measured in 
the lower resolution image (77 mJy), thus indicating the 
presence of a diffuse structure of low brightness that, ac- 
cording to Feretti et al. I (|2006f ) we associate to the relic 
structure. Because of the low flux density of discrete sources 
(6.9 and 4.2 mJy, respectively), the flux density and size of 
source C are only marginally affected by the presence of 
these possibly unrelated sources. This source could belong 
to the subclass of small-size relics like those in A13, A85, 
A133, and A4038 (Giovannini and Feretti 2004, Feretti and 
Giovannini 2007) located near the central brightest cluster 
galaxy. A proper study of this central region is prevented by 
the presence of nearby discrete sources and the insufflcent 
angular resolution of the X-ray data. 



4.2. Pressure in relic regions 

The Mach number obtained in the previous subsection can 
be used to estimate a theoretical pressure ratio P2/P1 ~ 
3.0. This is consistent with the observations: indeed in the 
direction relic A the pressure ratio is 2.9, of relic B it is 2.3, 
toward the N direction it is 2.8, and between A and B it is 
2.7. In the other directions the pressure ratio is much lower: 
to the W it is 1.6, and to the E it is 2.0 (see Tableg]). These 
results conflrm that the shock region is only in the relic 
direction. Note that to calculate the pressure, we used the 
observed electronic density in each direc tion obtained fr om 
/?- model flts of the surface brightness (^Neumann II2005I ). 



We computed the minimum nonthermal energy density 
in the relic sources from the radio data under equipartition 
conditions, with standard assumptions. We obtained about 
1.6 xlO~^^ erg cm~^ for relics A and B and 1.2x10"^^ erg 
cm~^ for C. The nonthermal pressure is 9.9 xlO~^^ erg 
cm~^ and 7.4 xlO~^^ erg cm~^, respectively. To estimate 
the minimum nonthermal energy density in the source C, we 
used the total flux density at low resolution, after subtract- 
ing components 7a and 7c (see Sect. 4.1). If we compare 
these results with the thermal pressure estimated in the 
relics region, we note that the minimum nonthermal pres- 
sure is about a factor 10 lower (but still in agreement with 
present uncertainties) than the thermal pressure. (From X- 
ray analysis we obtained 3.5x10"^^ erg cm~^ in the shock 
near relics A and B, and ~ 10~^^ erg cm~^ in the relic re- 
gion.) This result conflrms our general analysis of the gas 
pressure in this region. The higher value for relic C agrees 
with locat at the cluster center, and not simply projected 
onto it. The equipartition magnetic fleld is 0.9 /iG in relics 
A and B, and 3.6 /iG in source C. 



Table 4. Gas parameters in several regions, near the A, B 
relic location, inside and outside the probable shock. 





T 


rie 


P 




region 


(keV) 


(10-3cm-3) 


(10~^^dyn cm" 


-') 


to A, in shock 


5.7 ± 1.0 


0.40di 0.01 


3.6 




to A, out shock 


3.3 ±0.5 


0.24 zb 0.01 


1.2 




to B, in shock 


5.9 =b 1.0 


0.33 ± 0.01 


3.2 




to B, out shock 


4.7 =b 1.1 


0.19=b 0.01 


1.4 




N, in shock 


5.7 =b 0.6 


0.36=b 0.01 


3.3 




N, out shock 


3.5 =b 0.5 


0.22zb 0.01 


1.2 




W, in shock 


3.8 =b 0.5 


0.39=b 0.01 


2.4 




W, out shock 


3.8 ± 0.5 


0.24zb 0.01 


1.5 




E, in shock 


3.7 ± 0.4 


0.45=b 0.01 


2.6 




E, out shock 


3.3 ± 0.6 


0.26± 0.01 


1.3 




relic A, 9' 


6.2 ± 2.2 


0.14± 0.01 


1.4 




relic B, 9' 


3.7 ± 1.2 


O.Mzb 0.01 


0.9 




between AB, in 


5.7 =b 0.6 


0.37=b 0.01 


3.4 




between AB, out 


3.5 ± 0.5 


0.22± 0.01 


1.2 





5. Cluster dynamical state 

A548b is a cluster with complex dynamical state. 
Observational evidence suggests the presence of a shock 
structure before the radio relics A and B. To explain this 
result, we need to understand the cluster's dynamical his- 
tory better since radio relics would generally be expected 
to coincide with the shock region. 

We therefore performed a study of the geometry of the 
merger from optical analysis and from simulations. In this 
study, we need to take into account the projection effects. 

5.1. A548b galaxies velocities and environment 

To estimate possible projection effects, we retrieved op- 
tical information on the galaxies of A548 from the NED 
database. Among 606 galaxies found in NED near the 
A548b cluster, 371 have a velocity measurement. Among 
these galaxies, 193 lie at a projected distance less than 
R200 ?^1.5Mpc (i.e. r<30.50 from the A548b cluster X- 
ray center. In Fig. [9l we present the histogram of the 
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Fig. 9. Histogram of the NED galaxies with known veloc- 
ity at a projected distance less than R200 from the X-ray 
center. The light green line is the Gaussian fit. A large ay 
= 1300 kms~^ is necessary to fit a single Gaussian to the 
data. If one tries with two Gaussians (red continuum and 
dashed lines), the ay of the major component drops to ~ 
925 kms~^ and the two velocity centers found are very close 
to the velocity of A548-2 and A548a 



velocities of these galaxies. First, we obtained a fit to 
this histogram using a single Gaussian. The fit result im- 
plies a large <jy value of nearly 1300kms~^. Using this 
value, the expected cluster average temperature should be 
kT ^ 8 keV more than tw ice the observed value in the 
relation kT = {^^)^'^^ (|Xue et al. I I2OOOD . Therefore, 
even if statistics are relatively poor, we have tried to fit 
the galaxy velocity histogram with two Gaussians. In this 
case, we found the major component with ay 0^ 925 
±135kms~^ and z = 0.0432, whereas the second compo- 
nent has cry c^ 664±200kms-^ and z = 0.038. 

Remembering that A548b (05h45m27.9s -25d56ml9s V 
= 12708 kms~^ z = 0.042389) belong to a structure contain- 
ing A548-2 (05h42ml8.0s -26d05m00s V = 12861 kms"^ z 
= 0.0429), probably in the same plane of the sky and A548a 
(05h48m38.9s -25d28ml2s V = 11834 kms"^ z = 0.039474), 
some 1000 km s~^ less than the two others, the velocities 
found in the fit with two Gaussians, compare quite well 
with these redshifts (see Fig. [9j). 

We then divided the galaxy data set in two groups at a 
redshift limit of 0.04. We have 39 galaxies with z < 0.04 and 
96 galaxies with z > 0.042. In an Ra-Dec plane centered on 
the X-ray center, we plotted (Fig.fTOb) these two velocities 
groups as blue and red circles. One can see that the "red" 
group gently fills the whole space, while the "blue" one 
seems to be more distributed along the NE - SW direction. 

We now try to estimate the mass ratio of the merging 
units. On one hand, if we split the histogram in two parts at 
z :^0. 04-0. 042 and consider identical galaxy mass distribu- 
tions for the "blue" and "red" samples, then the number ra- 
tio is a first approximation of the mass ratio. In this way, we 
find :^2.46. On the other hand, if we use the fit results, the 
ratio between the areas of each component gives an estimate 
of 0^ 2.70±0.94 or, if we try with ay, and M1/M2 0^ (^)^ 
we find a mass ratio :^ 1.9±1.3. If we remember the poorf 
statistics, these three estimates are fairly coherent even if 
not strongly constraining. 



We can thus explain the apparently large ay and the 
spatial distribution of the "red" galaxies if we consider the 
following merging scheme: A548b would be the merger of 
two units, the major component being A548b, and a smaller 
one from the direction of A548a. Notice that, in this scheme, 
the collision axis forms a relatively small angle with respect 
to the line of sight. 

Different situations are possible as a merger between 
A548b and A548a, from the NE or clump arriving from 
A548-2 to A548b but it is very probably so that the merger 
collision is nearly perpendicular to the plane of the sky. 

In conclusion, from this optical analysis, we propose 
that we are observing a merger between two clusters with 
c^ 1:2-3 mass ratio (estimated from the velocity histogram 
of galaxies as shown in Fig. [9]). The merger's collision axis 
is nearly on the line of sight. The smaller cluster probably 
is located in the SW from the cluster emission center and 
we are observing the merger just after the core collapse. 
It should be noted that we detected no galaxy associated 
with the second component ("blue" circles), behind the ra- 
dio relics A and B, and in the cluster center (see Fig. [TQ|) . 

5.2. Comparison with simulations 

To understand the dynamical history of A548b and to ex- 
plain the connection between X-ray, radio, and optical data, 
we found it useful to compare observational data to the the- 
oretical predictions of a high-resolution cosmological simu- 
lation. For that purpose, we considered the 6th most mas- 
sive halo (hence the name Cluster6) in the z = snap- 
shot of a dark matter only simulation of a ACDM universe 
with a periodic box size of 80 h~^Mpc. We then resimu- 
lated the same box this time using the N body and hy- 
drodynamics AMR code RAMSES (Teyssier 2002), with a 
much higher resolution in a region of 20 h~^Mpc around 
the cluster, and with a much lower resolution outside this 
radius, just enough to get the correct large-scale tidal field. 
This now rather standard simulation technique is usually 
called "zoom simulation" or "re-simulation" , and it is used 
in cosmology to focus computational resources on one spe- 
cific object. Cluster6 is very interesting for this paper, 
since its final Virial temperature is close to 3.3 keV and 
R200 ~1.42Mpc, in striking agreement with A528b. Last 
but not least, Cluster6 experienced a major merger around 
z :^ 0.7, and at this redshift Cluster6 should be close to the 
assumed dynamical state of A548b. The simulation was ini- 
tialized using a 512^ Cartesian grid to cover the whole box, 
and then derefined outside the central 20 h~^ sphere. We 
then progressively and dynamically increased the level of 
refinement only in the high-resolution region^ using a quasi- 
Lagrangian refinement strategy, until we reached 7 addi- 
tional levels of refinement. Our finest cell size is therefore 
close to 1 /i~^kpc, comfortabily smaller than the observa- 
tional resolution. The number of AMR cells in the Virial 
radius is almost 1 x 10^, while the number of dark matter 
particles is slightly above 7 x 10^. 

We performed the comparison between simulation re- 
sults and X-ray data for gas density, gas temperature and, 
Mach number. These parameters are important for tracing 
the connection with radio data and particle acceleration. 

We chose the different outputs obtained from simula- 
tion around the maximum core collapse. We focused on the 
epoch when two clumps merge along a large-scale filament 
and consider different projection axes. At the same time. 
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Fig. 10. a) X-ray from ROSAT data and distribution of galaxies, the positions of the galaxies around the X-ray center, 
b) X-ray, radio data and "blue" circle correspond to galaxies having z < 0.04. c)X-ray, radio data, and "red" circles 
correspond to galaxies with z > 0.042. 






Fig. 11. Simulation data. We treated outputs just after 
maximum core collapse, in different projections. The figure 
shows the same output for Y projection and for Z projec- 
tion, in particular Y gas density map, Y gas temperature 
map with color scale, Z gas density map, Z gas tempera- 
ture map with color scale. We observed the shock in 210 
kpc with ratio T2/T1 = 2 in Y projection. 



we studied the effect of projection on the gas distribution. 
We examined the different projection X, Y, Z of the merger 
to have a complete scenario of all the merger situations. We 
obtained different results for the gas density, gas tempera- 
ture map with a the same outputs, from the same cluster at 
the same time. After our analysis of different outputs from 
simulation and comparison with X-ray data of A548b, we 
selected the output results just after the maximum core 
collapse. In this case simulated results for the temperature 
ratio and distance to the shock agree with observational 
results of A548b. 

Figure [11] shows the temperature and density maps ob- 
tained from the Y and Z projections. The results are very 
different, and we select the Y projection, because the den- 
sity map has the same distribution of the A548b emissivity 
map. Moreover, the temperature map obtained from the 



simulation has a temperature increase around the cluster 
center, similar to the observed one. The maximum of the 
temperature is in the annulus at 210 kpc, with values be- 
tween 4.5 and 6 keV. The same temperature of A548b (3.3 
keV) was obtained at the cluster center. Two peaks in tem- 
perature (6 keV) are present in the two regions around the 
cluster center. Note that the position of the peaks depends 
on the line-of-sight direction. After the optical analysis, we 
derived that A548b is in the state of merging, with 1:2-3 
mass ratio of the units, and the merger collision direction 
nearly perpendicular to the plane of the sky. 

To compare our results with radio data, we calcu- 
lated the Mach number map for different projections of 
chosen outputs, using the projected velocity distribution 
and the projected speed of sound, in particular (|Vazzal 
iBrunetti fc Gheller 2008h : 



Av 



3 1 



m2 



(3) 



Figure [12] shows the obtained Mach number maps for Y, Z 
projection, just after maximum core collapse. In Y projec- 
tion we observed a circular shock after temperature shock 
in the radius 250kpc, where the relativistic particles can be 
accelerated. 




Fig. 12. Obtained Mach number from simulation at the 
same output, just after maximum core collapse for the Y 
and Z projection. 



6. Discussion 

To explain the observational properties of A548b we need 
to understand the cluster dynamical history, in particular 
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Fig. 13. Probable geometry and merger collision scenario. 
Red regions indicate the shocks, green regions indicate the 
relics. The arrow indicates the merger direction. 

the merger collision scenario. After our analysis we found 
the following main facts: 

— From the X-ray analysis, we observed the shock in the 
ICM temperature and emission just before the relic's 
location, and the increase in the temperature and highly 
perturbed emissivity profiles in all cluster sectors. 

— From the optical analysis, we found that we very likely 
observed a merger of two units with 1:2-3 mass ratio, 
with the smaller unit located in the SW. 

— From our simulations and from A548b environments, we 
found that the merger collision should be nearly perpen- 
dicular to the plane of the sky. 

To reconstruct the geometry of the merger, we com- 
pared the observations with simulations. We propose the 
following merging scheme: A548b would be the merger of 
two units, with a 1:2-3 mass ratio, a major component com- 
ing from A548b, and a smaller one coming from the NE 
direction ( i.e. from A548a). The merger collision is nearly 
perpendicular to the plane of the sky, the small clump ar- 
rived from the NE, merged with the large clump, and now 
the gas from the small clump is located in the SW region, 
with respect to the cluster emission center. 

In the scenario that radio relics should be recently trig- 
gered by a shock wave for an efficient particles acceleration 
mechanism, the projection effects could help us to reconcile 
the relic position with the apparent position of the temper- 
ature enhancement. The A548b complex scenario, i.e. the 
relic's location (on the same side with respect to the cluster 
center) and the shock just before the relic's location, can be 
explained by the merger geometry sketched in Fig. [131 The 
two radio relics A and B, located at projected distances of 
500 and 430 kpc from the cluster center, lie at the bound- 
ary of the X-ray brightness distribution, but this is due to 
strong projection effects of a merger along the line of sight. 

The two relics represent two enhanced regions belong- 
ing to a large-scale shock, originated at the cluster center 
and moving to the peripheral region. A second symmetri- 
cal shock is expected by simulations, but because of mass 
ratio and merging geometry, it should be present in a pe- 
ripheral region in SW direction (see Fig. 13) with a very 
low density. It cannot be present in our observation be- 
cause of the chosen pointing center, but we do not expect 
that it could also be seen in a pointed XMM observation. 
According to present models of the formation of radio relics, 
we can imagine that a bubble of relativistic particles is ac- 
celerated by the shock originated by the major merger, the 
local magnetic field is compressed and amplified, originat- 



ing the radio emission. The origin of relativistic particles is 
unknown. We note that relativistic particles should not be 
present in the whole shocked region, but only in the region 
where we see the extended radio emission. This suggests 
that these relativistic particles, possibly originating in past 
AGN activity, were confined in floating discrete bubble re- 
gions. 

In this scheme, we observe the shock before the relic 
location due to the projection effects, because the merger 
collision is nearly perpendicular to the plane of the sky. Due 
to the poor statistics, the temperature in the relic regions is 
poorly constrained, but we expect these regions also to be 
shocked regions. In this scenario, the relics are at the same 
distance from the cluster center as the detected hot region, 
although the projected distance is smaller. The origin of 
relativistic particles responsible for the relic radio emission 
is puzzling. Indeed, it seems that the relativistic particles 
are not present in the fully shocked region, but only where 
radio emission is detected. 

7. Conclusions 

The main results of our study of the galaxy cluster Abell 
548b from XMM-Newton observations can be summarized 
as follows: 

X—ray analysis. We used XMM-Newton data to obtain 
X-Ray images and temperature distribution of the thermal 
ICM in A548b. We found that the /3-model is not a good fit 
of the data for the whole radial surface brightness profile in 
different sectors towards the N, S, W and E. Moreover we 
observed significant perturbations at 1.5', 5', and 7' from 
the cluster center in all the surface brightness profiles we 
obtained in different sectors and directions. The most sig- 
nificant perturbation is present towards relics A and B. 
Perturbations are also present in the residual image after 
subtraction of a 2D /? model, and here the most relevant 
perturbations are in the region towards the relics. 

We studied the dynamics of A548b from a detailed spec- 
tral analysis. The central cluster temperature estimated in 
the inner 3' is 3.4± 0.1 keV. The spectral study conffims 
that this cluster is not relaxed. We obtained a radial tem- 
perature profile, which shows a clear temperature increase 
at 5' from the cluster center. We extracted spectra in the 
regions between the cluster center and the relics A and B, 
and found that the most significant increase in temperature 
is present in the N (in direction of the two relics), but defi- 
nitely before the relic location. We studied the temperature 
distribution in sectors, in boxes, and with the X-ray wavelet 
spectral mapping algorithm. In all cases we obtained simi- 
lar results: a temperature increase in the annulus at 5' from 
the cluster center. The most significant increase in the tem- 
perature (from 3.5 keV to 6 keV) was observed towards the 
relic's direction. The temperature distribution shows a sim- 
ilar asymmetry to the brightness distribution. on. 

Optical analysis. From an optical analysis we observed 
a large galaxy velocity dispersion of cry :^ 1300kms~^, 
which is likely to be derived from two groups of galax- 
ies with different velocities and velocity dispersions: "red"- 
av ^ 935kms-^ and z = 0.0432; "blue"- ay ^ 664kms-^ 
and z = 0.038. Assuming that the "red" galaxies are associ- 
ated with the cluster A548b, the "blue" galaxies are associ- 
ated with a merging clump, with a 1:2-3 mass ratio. We can 
explain the apparent large <jy and the spatial distribution 
of the "blue" galaxies if we consider the following merging 
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scheme. Since A548b is in a rich environment with two other 
clusters A548-2 and A548a, it would be the merger of two 
units, a major component and a smaller one coming from 
the NE (from A548a). From the optical data we know, that 
A548a, A548b, and A548-2 are located at redshift 0.0394, 
0.0423, and 0.0429, respectively. In this scheme, the colli- 
sion axis forms a relatively large angle with the plane of the 
sky. 

Simulation analysis. We used Cluster 6, obtained from 
the hydro N-bo dy cosmology simulation, code RAMSES 
(|Tevssier Il2002f ). to draw the dynamical history of A548b. 
The mean temperature obtained from simulations is 3.3 
keV, consistent with the temperature observed in A548b. 
We compared simulations and X-ray data for density, tem- 
perature, and Mach number maps at epochs around the 
maximum core collapse, in particular just after the maxi- 
mum core collapse. We also studied the influence of projec- 
tion effects on the observations. After this simulation anal- 
ysis, we selected the output where density and temperature 
simulated maps were very similar to images obtained from 
our observation. The better agreeement was found under 
Y projection, i.e. with the merger axis nearly perpendicu- 
lar to the plane of the sky, and at the epoch just after the 
first maximum core collapse. To compare our results with 
the radio data, we also derived the Mach number map from 
the previous simulations. 

X — ray and radio connection. A548b is a cluster in a 
complex dynamical state. From the current data, it is not 
possible to establish if the central relic C is associated with 
any shock or X-ray feature. Higher angular resolution data 
are needed to analyze this cluster region. We observed an 
increase in the temperature at 5' -7' from the cluster center, 
i.e. before the radio relics A and B. This is interpreted as 
being caused by a shock with Mach number ~1.5. We can 
explain the relative location of the relics with respect to the 
shocks by strong projection effects due to a merger occur- 
ring nearly perpendicularly to the plane of the sky (see Fig. 
[T3|) . This scenario agrees with the su^^estio n that relics A 
and B could be part of a common structure (|Feretti et al. I 
I2006D . 
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